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Abstract: Urumqi Glacier No. 1 is a representative glacier in the inland areas of Central Asia and is the 
only Chinese reference glacier in the World Glacier Monitoring Service. In this study, we explored 
multi-decadal variations in the flow velocity of the glacier and the influencing factors based on continuous 
field observations and path coefficient analysis. Results show that the glacier flow velocity decreased from 
5.5 m/a in 1980/1981 to 3.3 m/a in 2010/2011. The annual variation in the direction of glacier flow 
velocity in the western branch and eastern branch was less than 1°—3°, and the change of glacier flow 
velocity in the western branch was more dramatic than that in the eastern branch. Glacier flow velocity 
was influenced by glacier morphology (including glacier area, glacier length, and ice thickness), glacier 
mass balance and local climate conditions (air temperature and precipitation), the glacier morphology 
being the leading factor. The long-term flow velocity data set of Urumqi Glacier No. 1 contributes to a 
better understanding of glacier dynamics within the context of climatic warming. 
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1 Introduction 


Glaciers and ice sheets are sensitive to climate change and preserve abundant water resources 
(Oerlemans and Fortuin, 1992; Haeberli et al., 2002; Oerlemans, 2005; Meier et al., 2007; Li et al., 
2012; IPCC, 2013). Global warming has driven a significant global reduction in the mass balance 
of most glaciers during the second half of the 20" century, resulting in a general retreating trend. 
This retreat was intensified at the end of the last century (Dyurgerov and Meier, 2000; Khromova 
et al., 2003; Paul et al., 2004; Barry, 2006; Li et al., 2010; Wang et al., 2013). Glacier flow 
velocity refers to the rate of downhill motion under the influence of gravity, and this motion 
process includes sliding on the glacier bedrock and the internal deformation of the glacier (Singh 
and Singh, 2001). Flow velocity is a basic parameter of glacier and is one of the main characters 
that differentiate glaciers from other natural ice bodies (Huang and Sun, 1982; Kääb, 2005). 
Glacier flow velocity varies at spatial and temporal scales. Monitoring of glacier flow velocity is 
a key to understanding the dynamics of glaciers (Zhou et al., 2009; Cuffey and Paterson, 2010). 
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Glacier flow velocity is a main factor controlling the glacier dynamics and depends on the 
dimensional parameters and physical properties of the glacier, such as ice thickness, surface and 
bed slopes, and ice temperature (Hodge, 1974; Mair et al., 2001). Under global warming, glacier 
loss can lead to changes in these parameters (or properties) and may then result in flow velocity 
change. 

In previous studies, various methods have been adopted to study glacier flow velocity, 
including field measurements (e.g., Manson et al., 2000; Hubbard and Glasser, 2005; Sugiyama et 
al., 2013) and remote sensing technologies (Strozzi et al., 2002; Nishimura et al., 2013; Schubert 
et al., 2013). Generally speaking, repeated field measurements of ablation stakes are the first and 
best choice for studying glacier flow velocity with high spatial and temporal accuracy (e.g., 
Manson et al., 2000; Hubbard and Glasser, 2005; Sugiyama et al., 2013), although repeated field 
measurements are only applicable to few easy-to-access glaciers and are quite labor and time 
consuming. 

Urumqi Glacier No. 1 is the only Chinese reference glacier in the World Glacier Monitoring 
Service (WGMS) and regarded as a representative for glaciers in the inland areas of Central Asia 
(Zemp et al., 2009). Flow velocity of Urumqi Glacier No. 1 has traditionally been measured by 
placing networks of stakes in the glacier surface and surveying their displacement at regular time 
intervals. Glacier flow velocity measurement was once taken in September 1959 and in 
August-September 1973. Since 1980, glacier flow velocity was systematically and persistently 
surveyed every year by various instruments, e.g., theodolite, total station and RTK-GPS (real time 
kinematic-global positioning system) (Tianshan Glaciological Station, 2011). In recent years, field 
surveys of stakes become easier and more accurate via the introduction of global positioning 
systems (GPS) (Zhou et al., 2009; Wang et al., 2014a). The monthly observation-based seasonal 
changes of glacier flow velocity of Urumqi Glacier No. 1 in 2006—2008 were previously reported 
(Zhou et al., 2009). However, long-term variations of glacier flow velocity have not yet been 
reported. 

The aims of this study are: (1) to analyze the spatial distributions of glacier flow velocity of 
Urumqi Glacier No. 1 in 1980/1981, 1990/1991, 2000/2001, and 2010/2011 with the time interval 
of 10 years, and (2) to present the long-term changes in glacier flow velocity from 1981 to 2011. 
The influencing factors are also discussed. 


2 Study area 


The Tianshan Mountains range is one of the major mountian ranges in Central Asia, with peaks 
rising from ~4000 to ~6000 m a.s.l. Urumqi Glacier No. 1 (43°06'30’N, 86°48'30"E; Fig. 1), a 
continental valley glacier, is located at the headwaters of Urumqi River in the Tianshan 
Mountains, Northwest China. The glacier covers northeast-facing slopes on the flanks of the 
Tianger Peak II with an altitude of 4484 m a.s.l. Both glacier accumulation and ablation occur in 
summer. The westerlies jet steam prevails across the area throughout the year. 

The area of Urumqi Glacier No.1 was reduced by 0.3 km? (15.7%) during 1962-2009. The 
terminus retreated with an average rate of 4.5 m/a from 1959 to 1993, and subsequently the 
glacier was separated into two independent branches (i.e., eastern branch and western branch) in 
1993. During 1994—2008, the terminus retreat rate was 6.0 and 3.5 m/a for the eastern branch and 
western branch, respectively (Li et al., 2011; Xu et al., 2011). According to the measurement 
results of B-1 radar-sounding system designed by Lanzhou Institute of Glaciology and 
Geocryology, Chinese Academy of Sciences (current name: Northwest Institute of 
Eco-Environment and Resources, Chinese Academy of Sciences) and ground penetrating radar, 
the average ice thickness was 55.1, 51.5, and 48.4 m for 1981, 2001, and 2006, respectively, and 
the bedrock elevation decreased towards the glacier terminus, exhibiting a U-shaped subglacial 
valley (Wang et al., 2014b). Li et al. (2012) developed an extended “perfect-plasticity” method to 
estimate the ice thickness along the glacier flow line of Urumqi Glacier No. 1. Comparision 
results of DEMs (Digital Elevation Models) from different periods showed that Urumqi Glacier 
No. 1 has experienced an accelerated thinning and a volume loss since the 1960s (Wang et al., 
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2014b). The accelerated mass loss was also observed from the annual observation of glacier mass 
balance since 1959 (Liu et al., 1997; Han et al., 2006; Jing et al., 2006; Li et al., 2011; Zhang et 
al., 2014). In addition, borehole drilling results showed that the ice temperature has increased by 
0.9°C (increasing rate of 0.06°C/a) during 1986-2001 and by 0.4°C (increasing rate of 0.08°C/a) 
during 2001—2006 (Li et al., 2011), further leading to the acceleration of glacier shrinkage. 
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Fig. 1 (a) Overview of Urumqi Glacier No. 1 and the surrounding area at the headwaters of Urumqi River, 
Tianshan Mountains, Northwest China, and (b) the surface elevation contours of Urumqi Glacier No. 1 in 2009 
and glacier boundaries in different years indicated by various color lines. It should be noted that a total of 7 
glaciers are located at the headwaters of Urumqi River, marked by No. 1-7, respectively. 


3 Data and methods 


3.1 Glacier flow velocity 


Glacier flow velocity was measured by surveying stakes installed in the glacier. A 42-stake 
network was set up across the entire glacier (Fig. 1b). These 42 stakes were evenly distributed in 
different elevations, covering nearly the whole glacier. It should be noted that only few stakes 
were distributed at the headwalls of the glacier due to inaccessible steep slopes. The velocity 
measurement stakes can also be used as mass balance measurement stakes. The total number of 
stakes varied from year to year, but it was never smaller than half of the recent number (Tianshan 
Glaciological Station, 2011). 

As mentioned above, sporadid glacier flow velocity measurements were taken in September 
1959 and also in August-September 1973 (Tianshan Glaciological Station, 2011). Since 1980, the 
flow velocity of Urumqi Glacier No. 1 was continuously surveyed in late August or early 
September every year. The measurements in different years were performed via various survey 
methods. Specifically, measurements of glacier flow velocity in 1959, 1973 and 1980-2001 were 
accomplished via repeated theodolite measurements from observation stations on the bedrock 
ridge near the northern edge of the glacier terminus. During 2002—2006, glacier flow velocity was 
measured using the total station. Since 2007, the RTK-GPS, a widely-used method for velocity 
measurements in recent years (Sunil et al., 2007; Rajner, 2010; Yang et al., 2014), was introduced 
to survey flow velocity. Our survey was undertaken using a Unistrong E650 GPS. Receivers were 
placed at a fixed base position that was the same as the previous year, and roving receivers were 
used to concurrently survey the stakes. The repeat measurements can be used to construct a flow 
velocity field. A RTK differential mode was simultaneously adopted, which resulted in a survey 
error of ~0.10—0.30 m for the geodetic-quality GPS receivers (Rivera et al., 2005). 

Furthermore, other various errors may also be introduced in the measurements of glacier flow 
velocity. First, some stakes fell down or moved from their original positions due to the 
complicated morphology, rapid glacier flow velocity, or severe glacier melting. Therefore, those 
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stakes need to be reinstalled to the original positions. However, sometimes the needed 
reinstallation is impossible under difficult situations (e.g., movement into crevassed areas or 
glacier surface streams). Second, the measurement duration must be sufficient to assure the 
survey accuracy, otherwise, errors may occur. For example, the GPS measurement for each stake 
requires at least 5 min. 

Glacier flow velocity was obtained by calculating the displacement of the stakes per unit time. 
The coordinate system was independent. The displacement of a stake was determined by taking 
the difference between each component of the coordinate. Then by dividing the displacement by 
the unit time, the glacier flow velocity in each direction per unit time was acquired. The 
magnitude of glacier flow velocity (Uxy, in unit of m/a) and the direction of glacier flow velocity 
(a, in unit of °) were determined using the following equations: 


U „=U? +U, (1) 


a soj- (2) 
where Ux and U; are the flow velocity components parallel to the x and y axes (m/a), respectively; 
and a is the direction of glacier flow velocity (°), which is the degree relative to north. Generally, 
the accuracy of U,, is less than 10% of the error in the input data and the accuracy of a is no more 
than +2° (Tianshan Glaciological Station, 2011). 

To examine the variations in flow velocity of Urumqi Glacier No. 1 from 1981 to 2011, we 
analyzed the spatial distributions of glacier flow velocity in 1980/1981, 1990/1991, 2000/2001, 
and 2010/2011 and the changes of glacier flow velocity with the time interval of 10 years. In 
addition, we chose the topographic maps containing Urumqi Glacier No. 1 made in 1981 (1:5000; 
aerial photograph), 1986 (1:5000; stereophotogrammetric survey), 2001 (1:5000; theodolite), and 
2009 (1:5000; RTK-GPS) to determine the glacier boundary in these years. More details about 
these topographic maps can be found in Wang et al. (2014a). 


3.2 Factor analyses 


Glacier flow velocity may be influenced by climatic condition, glacial morphology, glacial mass 
balance, subglacial topography, subglacial water pressure, and internal stress (McSaveney and 
Gage, 1968; Iken and Truffer, 1997; Herman et al., 2011; Nishimura et al., 2013; Sugiyama et al., 
2013). In this study, the correlation coefficients between various variables were determined by 
statistical analysis in the software Microsoft Office Excel. The P-value was used to test the 
statistical significance and the difference was considered to be significant at P<0.05 level. 
Furthermore, we adopted the path coefficient analysis to determine the direct and indirect 
influences of causal variables on glacier flow velocity. Developed by Wright (1921), the path 
coefficient analysis is a method that quantifies the contribution of causal variables to a targeted 
influence variable. This method was originally used in the context of genetics (Dewey and Lu, 
1959). It has since been widely applied to many other fields (e.g., Board et al., 1997; Machikowa 
and Laosuwan, 2011; Akintunde, 2012). The method can be expressed using the following 
formula by assuming independent variables x,, x2, x3, ..., x, and dependent variable y (Jing et al., 
2011): 


Pyx, +1 x, Pyx, be thy, Px, = ry 


Pex Pyx, + Pyx, Hers Tex, Pyx, = Foy 


; (3) 
Pagpe, P yxy Y Txpxs Pyrs TUE Py, gy 


where, r... 
itj 


is the correlation coefficient between x, and x Nay is the correlation coefficient 
between x; and y; p,, is the direct path coefficient of x; and y; and Tx, Pye is the indirect path 


coefficient of x; and y via x. 
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The direct path coefficient ( p,,.) can be expressed as Equation 4: 
Py, =0,6,,/06, ((=1,2,3,....7), (4) 
where, b; is the partial regression coefficient of y and x; and Oy and ô, are the standard 


deviations of x; and y, respectively. The path coefficient of residual terms can be expressed as 


Equation 5 and the large value indicates the deviation is relatively large or some important factors 
are not considered. 


Pyx, = EEN +1, yPyx, +7 tT, yPyx,,)+ (5) 


In our study, glacier flow velocity was chosen as the dependent variable, while glacier area, 
glacier length, glacier mass balance, ice thickness, air temperature and precipitation were selected 
as the independent variables. The data of glacier area and length have been published in our 
previous study (Wang et al., 2014a). It should be noted that ice thickness of Urumqi Glacier No. 1 
was measured in 1981, 2001, and 2006 (Wang et al., 2014b), and the ice thickness used for other 
years (1981, 1986, 1992, 1994, 2000, 2001, 2006, and 2009) was calculated by the empirical 
formula for uniformity (Liu et al., 2003). Glacier mass balance was observed by the stake/snow 
pit method since 1959. Air temperature and precipitation data were obtained from the Daxigou 
Meteorological Station (43°06'N, 86°50’E; 3539 m a.s.l.), which is 3 km southeast of Urumqi 
Glacier No. 1. A continuous time series of air temperature and precipitation data are available 
from the station, which has been operated by the Xinjiang Uygur Autonomous Region 
Meteorological Bureau. 


4 Results and discussion 


4.1 Spatial and temporal distribution of glacier flow velocity 


Glacier flow velocity values in 1980/1981, 1990/1991, 2000/2001, and 2010/2011 are shown in 
Figure 2. For the years 1980/1981, 1990/1991, 2000/2001, and 2010/2011, the values of average 
glacier flow velocity were 5.5, 4.6, 3.8, and 3.3 m/a, respectively; and the values of maximum 
glacier flow velocity were 10.6, 8.3, 6.1, and 4.6 m/a, respectively (Fig. 2a). The decreasing rate 
of average glacier flow velocity was 0.07 m/a from 1980/1981 to 2010/2011. Spatially, the flow 
velocity tended to decrease in the upper and lower parts of the glacier. Glacier flow velocity was 
largest near the main flow line and smallest near the margin. The bedrock of Urumqi Glacier No. 
1 shows a U-shaped pattern and is quite regular with altitudes decreasing towards the ice front 
(Wang et al., 2014b). The directions of glacier flow velocity in different years were all parallel 
with the valley direction. Generally, annual variations in the directions of glacier flow velocity 
were small, less than 1°—3°. 

Variations of glacier flow velocity for the western and eastern branches showed a consistent 
pattern from 1980/1981 to 2010/2011 (Figs. 2b and c). The average glacier flow velocity of the 
western branch was 1.4 times higher than that of the eastern branch. This may be due to the larger 
thickness and the steeper slope of the western branch (Wang et al., 2014b). The slope of the 
western branch is on average 8° higher than that of the eastern branch. The steeper slope can 
induce the ice to transfer from the accumulation zone to the ablation zone at a higher velocity. 
Assuming that ice deforms in simple shear fashion, glacier flow velocity (u,) can be given by the 


following equation (Cuffey and Paterson, 2010): 
2A : 
u, =u, +——(pgH sina)" H, (6) 
n+1 


where u, is the glacier flow velocity (m/a); u, is the basal velocity (m/a); p is the ice density (m/a); 
g is the gravitational acceleration (m/s?); H is the ice thickness (m); a is the down-glacier surface 
slope (°); and A is the rate factor and n is the exponent, both of them are parameters in Glen’s 
flow law. The second term of Equation 6 is proportional to H* and (sina)*, which states that the 
glacier flow velocity is highly dependent on the ice thickness and surface slope. 
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Fig. 2 Average and maximum flow velocities of the entire Urumqi Glacier No. 1 (a), the western branch (b), and 
the eastern branch (c) from 1980/1981 to 2010/2011 


Based on the ice thickness results of Urumqi Glacier No. 1 from ground penetrating radar 
survey (Wang et al., 2014b), we determined that locations with the maximum glacier flow 
velocity for the two branches were at the places where ice is the thickest (Fig. 3). The locations 
with the maximum glacier flow velocity for both branches moved backward (i.e., towards 
upstream) from 1980/1981 to 2010/2011. This is likely because the glacier thinned more in the 
lower part than in the upper part, resulting in a decrease of glacier flow velocity in the lower part. 

Urumqi Glacier No. 1 also showed seasonal variations in flow velocity (Zhou et al., 2009). 
Monthly observations of glacier flow velocity from 2006 to 2008 indicated that glacier flow 
velocity in summers was higher than that in winters. This may be due to the fact that the resulted 
melting water from high temperature in summers may have lubricated the ice-bed and speeded up 
the glacier flow. 

Compared with other glaciers in China, the flow velocity of Urumqi Glacier No. 1 is lower 
because of the smaller size. For examples, the average flow velocity of Qiyi Glacier in the central 
Qilian Mountains was 9.02 m/a from 1958 to 2007 (Jing et al., 2011). The maximum flow 
velocity reached up to 32.6 m/a in the western branch and 32.4 m/a in the eastern branch in 
Laohugou Glacier No. 12 (39°26'N, 96°32’E) in the western Qilian Mountains (Liu et al., 2010). 
In contrast, the maximum flow velocity of Yanglonghe Glacier No. 5 (39°13'N, 98°33’E) in the 
central Qilian Mountains was 8.8 m/a, and the maximum flow velocity of Chongce Ice Cap 
(35°14'N, 81°07’E) in the western Kunlun Mountains was 4.16 m/a. The maximum velocity of 
Puruogangri Ice Field (33°44'-34°04'N, 89°20'-89°50’E), the largest ice field in the northern 
Tibetan Plateau, was 3.67 m/a in 2000 (Jing et al., 2003). The thinning of Urumqi Glacier No. 1 
was remarkable during the period of 1962-2009 (Wang et al., 2014a), resulting in the decreasing 
of glacier flow velocity. During the period of 1959-2010, the cumulative glacier mass balance of 
Urumqi Glacier No. 1 was —14,883 mm water equivalent. And the equilibrium line altitude 
showed a significant ascending trend with a value of approximately 90 m lowering (Zhang et al., 
2014). 


4.2 Factor analysis 


Glacier flow velocity is mainly influenced by glacier slope, ice temperature, and ice thickness 
(Hodge, 1974; Mair et al., 2001). Generally speaking, glacier slope represents the morphological 
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Fig. 3 Glacier flow velocity vectors in (a) 1980/1981, (b) 1990/1991, (c) 2000/2001, and (d) 2010/2011. The 
glacier boundaries were determined using the topographic maps in 1981, 1986, 2001 and 2009, with the survey 
methods of aerial photograph, stereophotogrammetric survey, theodolite, and RTK-GPS, respectively. 


characteristics of glacier. Ice temperature is affected by air temperature to some extent. Moreover, 
ice thickness change is normally reflected by glacier mass balance and closely related to climate 
factors (air temperature and precipitation). Therefore, we analyzed the influences of glacier area, 
glacier length, glacier mass balance, ice thickness, air temperature and precipitation on glacier 
flow velocity using the path coefficient analysis in this study. 

Correlation coefficients among different variables are shown in Table 1. Glacier flow velocity 
was significantly correlated with glacier area (0.89), glacier length (0.92), glacier mass balance 
(0.90), ice thickness (0.89), and air temperature (0.80). Linear analyses of air temperature and 
precipitation data from Daxigou Meteorological Station (Fig. 4) showed that annual mean 
temperature increased by 0.28°C/10a and annual precipitation increased by 16.93 mm/10a from 
1959 to 2014. The increasing tendency was particularly obvious after 1980 with the rates of 
0.54°C/10a for annual mean temperature and 33.36 mm/10a for annual precipitation. Glacier flow 
velocity exhibited a significant negative correlation with air temperature (—0.80) and a 
non-significant negative correlation with precipitation (—0.26) (Table 1). It should be noted that 
because the ice thickness was calculated by the empirical formula, correlation coefficients 
between glacier flow velocity and glacier area and between glacier flow velocity and ice thickness 
were the same. Glacier mass balance showed a significant negative correlation with air 
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temperature (correlation coefficient of —0.89), which is similar to the result (correlation 
coefficient of —0.69) reported by Zhou et al. (2010). 


Table 1 Correlation coefficients between dependent and independent variables by statistical analysis 


GA GL IT GMB Ta Pre GFV 
GA 1.00 
GL 0.97” 1.00 
IT 0.99% 0.97" 1.00 
GMB 0.95" 0.98" 0.95” 1.00 
Ta -0.79% -0.84™ -0.79" -0.89" 1.00 
Pre —0.53 —0.48 —0.53 —0.40 0.04 1.00 
GFV 0.89" 0.92" 0.89" 0.90" -0.80" —0.26 1.00 


Note: GA, glacier area; GL, glacier length; IT, ice thickness; GMB, glacier mass balance; Ta, air temperature; Pre, precipitation; GFV, 
glacier flow velocity. ** indicates significant difference at P<0.01 level. Glacier area, glacier length, glacier mass balance, ice thickness, 
air temperature and precipitation are independent variables, while glacier flow velocity is the dependent variable. 
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Fig. 4 Changes of annual mean temperature and annual precipitation of Daxigou Meteorological Station from 
1959 to 2014 


Again, glacier flow velocity was significantly correlated with all the aforementioned variables 
and all the variables are independent (with an exception of precipitation; see Table 1). We thus 
used the path coefficient analysis to determine the direct and indirect influences of these factors 
on glacier flow velocity. The direct and indirect path coefficients between glacier flow velocity 
and those independent variables (glacier morphological factors, glacier mass balance, and climate 
factors) are shown in Table 2. The direct path coefficients of glacier area, glacier length, glacier 
mass balance, ice thickness, air temperature, and precipitation with glacier flow velocity were 
59.97, —1.44, 3.13, —59.89, 0.59, and 0.54, respectively. The path coefficient of residual terms was 
0.02, indicating that the main factors influencing glacier flow velocity have already been 
considered in this study. To sum up, the direct path coefficient of glacier morphological factors 
(glacier area, glacier length, and ice thickness) with glacier flow velocity was the highest, 
followed by glacier mass balance and climate factors (air temperature and precipitation) (Table 2), 
indicating that glacier morphology exerted the greater influence on glacier flow velocity than 
climatic conditions, being similar to the results from the Qiyi Glacier (39°14'13"N, 97°45'20"E) 
in the central Qilian Mountains, northern edge of the Qinghai-Tibetan Plateau (Jing et al., 2011). 
It was previously reported that glacier length and area play an important role in glacier flow 
velocity (Jing et al., 2011). For Urumqi Glacier No. 1, glacier flow velocity decreased with the 
reductions in glacier area and also with the retreats in terminus position. It should be particularly 
noted that more complex interactions may exist among these influencing factors and a thorough 
analysis is needed in future studies. 
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Table 2 Path coefficients between glacier flow velocity and various influencing factors 


; Correlation coefficient Direct path Indirect path coefficient 
Variable ; is 
with GFV coefficient GA GL IT GMB Ta Pre Total 
GA 0.89 59.97 -1.39 -59.89 2.96 -0.47 -0.29 -59.08 
GL 0.92 -1.44 58.10 —58.07 3.07 -0.50 -0.26 2.36 
IT 0.89 —59.89 59.96 -1.39 2.97 —0.47 -0.29 60.80 
GMB 0.90 3.13 56.77 -141 -56.84 —0.52 -0.22 -2.22 
Ta —0.80 0.59 —47.31 1.21 47.46 -2.77 0.02 -1.39 


Pre —0.26 0.54 -31.94 0.69 31.68 -1.26 0.02 —0.80 


5 Conclusions 


In this study, spatial distributions and long-term changes of glacier flow velocity in Urumqi 
Glacier No. 1, as well as the influencing factors, were analyzed. The average flow velocity in the 
western branch of Urumqi Glacier No. 1 was approximately 1.4 times faster than that in the 
eastern branch. During the study period (i.e., 1980/1981—2010/2011), the decreasing rate of 
average glacier flow velocity was 0.07 m/a and the decreasing rate may be directly or indirectly 
associated with the glacier retreating and thinning. The annual variation in the directions of 
glacier flow velocity was less than 1°—3°. Glacier flow velocity was primarily influenced by 
glacier morphology (glacier area, glacier length, and ice thickness), and the glacier mass balance 
and climate conditions (air temperature and precipitation) are also significant in influencing 
glacier flow velocity. For Urumqi Glacier No. 1, continuous field observations (e.g., flow velocity, 
mass balance, and terminus change) are needed in future to further our understanding of the 
glacier dynamics in the context of climatic warming. 
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